Spherical Dust Collapse in Higher Dimensions 
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We consider here the question if it is possible to recover cosmic censorship when a transition is 
made to higher dimensional spacetimes, by studying the spherically symmetric dust collapse in an 
arbitrary higher spacetime dimension. It is pointed out that if only black holes are to result as 
end state of a continual gravitational collapse, several conditions must be imposed on the collapsing 
configuration, some of which may appear to be restrictive, and we need to study carefully if these 
can be suitably motivated physically in a realistic collapse scenario. It would appear that in a 
generic higher dimensional dust collapse, both black holes and naked singularities would develop as 
end states as indicated by the results here. The mathematical approach developed here generalizes 
and unifies the earlier available results on higher dimensional dust collapse as we point out. Further, 
the dependence of black hole or naked singularity end states as collapse outcomes, on the nature 
of the initial data from which the collapse develops, is brought out explicitly and in a transparent 
manner as we show here. Our method also allows us to consider here in some detail the genericity 
and stability aspects related to the occurrence of naked singularities in gravitational collapse. 

PACS numbers: 04.20.Dw, 04.20.Cv, 04.70.Bw 



I. INTRODUCTION 

A considerable debate has continued in recent years on 
the validity or otherwise of the cosmic censorship conjec- 
ture (CCC) in black hole physics (see e.g. [1-7] for some 
recent reviews) . The reason for this interest is that CCC 
is fundamental to many aspects of theory and applica- 
tions of black hole physics, and the astrophysical impli- 
cations resulting from a continual gravitational collapse 
of a massive star which has exhausted its nuclear fuel. As 
of today, no theoretical proof or any satisfactory mathe- 
matical formulation of CCC is available, where as many 
collapse scenarios have been analyzed and worked out 
wherein the collapse end state is either a black hole (BH) 
or a naked singularity (NS), depending on the nature of 
the initial data from which the collapse evolves develop- 
ing from a regular initial state to the final super-dense 
state. This has important astrophysical significance for 
the reason that naked singularities may have observa- 
tional properties which could be radically different from 
those of a black hole. 

A note-worthy suggestion that has emerged towards 
a possible theoretical formulation of CCC is that any 
naked singularities resulting from matter models which 
may also develop singularities in special relativity, should 
not be regarded as physical (see e.g. [3]). Clearly, it will 
require a serious effort to cast such a possibility into a 
mathematical statement and a proof for CCC. Also, it 
may not be easy to discard completely all the matter 
fields such as dust, perfect fluids, and matter with vari- 
ous other reasonable equations of state, which have been 
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studied and used extensively in relativistic astrophysics 
for a long time. 

Another possibility that indeed appears worth explor- 
ing is we may actually be living in a higher dimen- 
sional spacetime universe. The recent developments in 
string theory and other field theories strongly indicate 
that gravity is possibly a higher dimensional interaction, 
which reduces to the general relativistic description at 
lower energies. Hence, while CCC may fail in the four- 
dimensional manifold of general relativity, it may well 
be restored due to the extra physical effects arising from 
our transition itself to a higher-dimensional spacetime 
continuum. Such considerations would inspire a study 
of gravitational collapse in higher dimensions. From 
such a perspective, many papers have reported results 
in recent years on spherically symmetric collapse of dust 
(where the pressures are identically taken to be zero) 
in higher dimensions [8-11]. The recent revival of in- 
terest in this problem is motivated to an extent by the 
Randall-Sundrum brane- world scenario [12]. Various au- 
thors considered specialized and different subcases of the 
general problem of dust collapse in higher spacetime di- 
mensions. For example, the marginally bound case in a 
general spacetime dimension was studied in [9], whereas 
the same case was studied in [10], however, with an addi- 
tional and physically motivated restriction on the initial 
density profiles, namely that the first derivative of the 
initial density distribution for the collapsing cloud must 
be vanishing. Also, the non-marginally bound case, how- 
ever, with the geometric assumption that spacetime is 
self-similar, was examined in [11] for a five-dimensional 
model. 

Such studies do provide us with an idea of what is 
possible in gravitational collapse as far as its end state 
spectrum is concerned. It is obvious from the work in 
this area so far that any possible proof of CCC must be 
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inspired by such additional physical inputs (one of these 
being a possible transition to higher dimension) into our 
current framework of thinking. Any such alternatives 
would be worth exploring due to the fundamental sig- 
nificance of CCC in black hole physics. The point is, if 
naked singularities did indeed develop in realistic gravita- 
tional collapse of massive objects, they may have proper- 
ties which would be rather different from those of black 
holes both theoretically as well as observationally, and 
a comparison of these two cases may prove to be quite 
interesting. 

From such a perspective, we investigate here the is- 
sue of BH/NS end states in a higher dimensional dust 
collapse in some detail, in order to bring out in a trans- 
parent manner the effect of dimensions on the final fate 
of the evolution of the matter cloud which collapses from 
a given regular initial data. A spherically symmetric col- 
lapse is considered in N > 4 dimensions and the matter 
form is chosen to be dust as this is a well-studied case. In 
fact, there have been suggestions in the past that at the 
vicinity of the collapse, the in-falling velocity of the mat- 
ter shells could be so high that the effects of pressures are 
negligible and hence dust may be a good approximation 
|13|. Apart from that, dust collapse is worth investiga- 
tion as it has continued to serve as a basic paradigm in 
black hole physics. 

As compared to some of the studies of dust collapse 
in higher dimensions mentioned above, we deal here 
within a general framework, i.e. there is no restriction 
adopted on dimensions, the spacetime is not assumed 
to be self-similar, which is a somewhat restrictive con- 
dition for collapse models, and we treat marginally and 
non-marginally bound cases together by means of a uni- 
fied treatment. Also, we make no assumptions on the 
nature of the initial data, i.e. the initial density and ve- 
locity profiles apart from their being regular, and these 
are taken to be suitably differentiable functions of the 
comoving radial coordinate at the origin. Our results 
thus generalize earlier results on the topic and provide a 
unified treatment. An additional advantage that is de- 
rived due to the treatment here is that the dependence 
of the BH/NS end states on the nature of initial data 
from which the collapse evolves, is brought out rather 
clearly. The methodology used also allows us to consider 
in some detail the genericity and stability aspects related 
to the occurrence of naked singularities in gravitational 
collapse. 

To predict the final state of collapse for a given ini- 
tial mass and velocity distribution, we study here the 
singularity curve resulting from the collapse of succes- 
sive matter shells, the tangent to which at the singu- 
larity is related to the radially outgoing null geodesic 
equation. Hence, by determining the sign of the tan- 
gent of the singularity curve at the central singularity 
we can integrate the null geodesies close to the singular- 
ity, and find whether a future directed radially outgoing 
null geodesic comes out from the shell-focusing central 
singularity. This allows us to determine the final fate of 



collapse in terms of either a black hole, or a naked sin- 
gularity. It turns out that for CCC to be true even in 
higher dimensions, one must impose a number of con- 
straints for the given system. The physical relevance of 
such assumptions is discussed here in some detail. 

The paper is organized as follows. In the next sec- 
tion we set up the Einstein equations, and the regularity 
conditions for collapse are specified. In Section III, con- 
tinually collapsing matter clouds are considered and it 
is shown how the initial data in terms of the initial dis- 
tributions of the density and velocity profiles determine 
the final fate of the collapse in terms of either a black 
hole or a naked singularity. The dependence of these 
BH/NS phases on the nature of initial data is brought 
out clearly here. In Section IV we discuss cosmic cen- 
sorship in higher dimensions from the perspective of our 
results, and the genericity and stability aspects related to 
naked singularity formation are also considered in some 
detail. The case of marginally bound collapse is discussed 
in some detail here together with some other cases of in- 
terest, relating this to the question of when it may be 
possible to restore the cosmic censorship for such models 
in terms of certain possible physical assumptions for the 
initial density profiles. Some conclusions and remarks are 
summarized in the final Section V. 



II. EINSTEIN EQUATIONS AND REGULARITY 
CONDITIONS 

To study the collapse of a spherical dust cloud, we 
choose a general spherically symmetric co-moving metric 
in N > A dimensions which has the form, 



ds 2 



-dt 2 + e 2 ^dr 2 + R 2 (t,r)dn 2 N _ 2 



where, 



N-2 

dfl 2 N _ 2 = ^ 

i=l 



(d6 i 



(1) 



(2) 



is the metric on (N — 2) sphere. The energy-momentum 
tensor of the dust has the form, 



P(t, 



t; = o ; t" 







(3) 



We take the matter field to satisfy the weak energy con- 
dition, i.e. the energy density measured by any local 
observer be non-negative, and so for any timelike vector 
V 1 . we must have, 



T ik V l V k > 



P > 



(4) 



In the case of a collapsing cloud, the cloud has a finite 
boundary, < r < rj,, outside which it is matched to a 
Schwarzschild exterior. The range of the coordinates for 
the metric is then < r < oo, and — oo < t < t s (r) where 
t s {r) corresponds to the singular epoch R — 0. 
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Using the above conditions, we can write the N- 
dimensional Einstein equations as [T^j. 

(N - 2) {^5^(1 + R 2 - R 2 e^') + 



1 , • 



R 



tyR + e-^ij}'R' - e'^R") \ = p (5) 



(N-2)\^(R'2e-^-R 2 -l) + H=0 (6) 



(N-2) 
R 



{2{R'-i,R')) 



= 



G\ (i = 2,--- ,N-2) = 



(7) 



(8) 



where Gf is the N-dimensional Einstein tensor. Inte- 
grating the first part of equation Q we get, 

R >2 



,2V 



1 + f(r) 



(9) 



where /(r) is an arbitrary function of co-ordinate r, 
and /(r) > — 1. Hence the generalized Tolman-Bondi- 
Lemaitre (TBL) metric in N-dimensions becomes, 



ds 2 = dt 2 



R 



-dr 2 - R 2 {t,r)dn 2 N _ 2 (10) 



1 + fir) 

Now substituting equation©, in (J5J and JSJ, we have, 

P, ^7^5 = (ID 



2R(. N ~ 2 )R' rl 2R( N ~VR 
where we define and get by solving eauation llljl . 

U = R (N - 3) (R 2 - f{r)) ; U = F{r) (12) 

Here F(r) is another arbitrary function of the comoving 
coordinate r. In spherically symmetric spacetimes F(r) is 
the mass function which describes the mass distribution 
of the dust cloud and f(r) is the energy function for the 
collapsing shells. Thus using equations (fTTft and (jT2"|) . we 
finally get the required equations of motion as, 



(N-2)F' _ 
2R( N - 2 )R' ~ 9 



^ 2 = ]|S)+/W (13) 



As seen from the above, once the mass function is spec- 
ified at the initial epoch the energy function f(r) fully 
specifies the velocity distribution of infalling shells. Also, 
the energy condition then implies F' > 0. It follows from 
the above that there is a spacetime singularity at R = 
and at R' = 0. The latter are called 'shell-crossing sin- 
gularities' which occur when successive shells of matter 



cross each other. These have not been considered gen- 
erally to be genuine spacetime singularities, and possible 
extensions of spacetime have been investigated through 
the same [l5|. On the other hand, the singularity at 
R = is where all matter shells collapse to a zero physical 
radius, and hence has been known as a 'shell-focusing sin- 
gularity'. The nature of this singularity has been inves- 
tigated extensively in four-dimensional spacetimes (see 
e.g. references in [4], and other reviews), and it is known 
for the case of spherical dust collapse that it can be both 
naked or covered depending on the nature of models one 
is considering. 

The collapsing matter cloud condition implies that 
R < 0. For dust clouds it follows form the equations of 
motion that once R < at the initial epoch from where 
the collapse commences, then at all epochs we have the 
same condition holding, and thus there is an endless col- 
lapse till the shell-focusing singularity R — is reached. 
In other words, there is no bounce possible in dust col- 
lapse once the collapse has initiated with R < 0. We use 
the scaling independence of the comoving coordinate r to 
write (see e.g. [4]), 



where, 



R(t, r) = rv(t, r) 



v(ti,r) = l ; «(«,(»■), r)=0 



v < 



(14) 



(15) 



This means we have scaled the coordinate r in such a 
way that at the initial epoch R = r and at the singu- 
larity, R — 0. It should be noted that we have R = 
both at the regular center r = of the cloud, and at the 
spacetime singularity where all matter shells collapse to 
a zero physical radius. The regular center is then distin- 
guished from the singularity by a suitable behaviour of 
the mass function F(r) so that the density remains finite 
and regular there at all times till the singular epoch. The 
introduction of the parameter v then allows us to distin- 
guish the spacetime singularity from the regular center, 
with v = 1 all through the initial epoch, including the 
center r — 0, which then decreases monotonically with 
time as collapse progresses to value v = at the singu- 
larity R = Q. 

From the equations of motion it is evident that to have 
a regular solution over all space at the initial epoch, the 
two free functions F(r) and f(r) must have the following 
forms, 



F(r) = r (JV4) jW(r); f(r) = r 2 b(r) 



(16) 



where M. (r) and b(r) are at least C 1 functions of r for r = 
0, and at least C 2 functions for r > 0. This is dictated by 
the condition that the density and energy distributions 
must be regular at the initial epoch and should not be 
blowing up. This is because if the mass function F did 
not go as power at least r*" 1 closer to the origin, then 
as implied from the equations of motion, the density will 
be singular at the origin r = as it will diverge there, 
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and that cannot be accepted as regular initial data for 
collapse. Similarly, equation (12) implies that f(r) is 
determined once the velocity profile is specified, and vice- 
versa, for a given initial density distribution. Since the 
center of the cloud is taken to be at rest in any spherically 
symmetric distribution, the leading term in the energy 
profile must be at least r or higher. Then again equation 
(12) implies the behaviour for f(r) as above. 



III. CONTINUAL COLLAPSE AND BH/NS 
END STATES 

We now consider the endless collapse of a dust cloud 
to a final shell-focusing singularity at R = 0, where all 
matter shells collapse to a zero physical radius. In par- 
ticular, we analyze specifically the nature of the central 
singularity at R = 0, r = in detail to determine when 
it will be covered by the event horizon, and when visible 
and causally connected to outside observers. If there are 
future directed families of non-spacelike curves, coming 
out from the singularity and reaching faraway observers, 
then the singularity will be naked. The absence of such 
families will give a covered case when the result is a black 
hole. 

In the following, a higher dimensional collapse as dis- 
cussed above is considered and we shall show how the 
initial density and energy distribution prescribed at the 
initial epoch from which the collapse commences, com- 
pletely determine if there will be families of non-spacelike 
trajectories coming out of the singularity. We shall show 
that for a generic situation, given an initial density dis- 
tribution for the collapse to develop, it is always possible 
to choose an energy profile so that the collapse of this 
density profile ends in a naked singularity, or one could 
also choose another class of energy distribution so that 
the same density distribution will end up collapsing so as 
to create a black hole. That is, given an initial density 
profile, the outcome in terms of either a black hole or 
a naked singularity as end states really depends on the 
class of the energy distribution chosen. The converse will 
also be seen to hold true, that is, given the initial energy 
function, one will be able to choose classes of density pro- 
files, subject to the the weak energy condition, so as to 
give rise to either of the BH/NS end states depending on 
the choice made. Our results reduce to the usual dust 
collapse models of general relativity when N = 4. 

With the regular initial conditions as above, equa- 
tion HI 3|) can be written as, 



v^v = -\ I M{r) + -y( Ar - 3 )6(r) 



(17) 



Here the negative sign implies that v < 0, that is, the 
matter cloud is collapsing. Integrating the above equa- 
tion with respect to v, we get, 



t(v,r) = 



Af-3 

v 2 dv 



v ^jM{r)+v( N -3)b(r) 



(18) 



We note that the co-ordinate r is to be treated as a con- 
stant in the above equation. Expanding t(v, r) around 
the center, we get, 

t(v, r) = t(v, 0) + rX(v) + r 2 ^ + + ■ • ■ (19) 

2 6 

where the function X(y) is given by, 

*(„) = -■ f'^(M 1+ ^%)^ 
2 /. (vWo + t)(~- 3 >& )* 

where, 

b = b{0) ; M o = M{0) ; h = b'(0) ; Mi = M'(0) 

(21) 

Thus, the time taken for the central shell to reach the 
singularity is given as 



ts — 



v 2 dv 



o y/Mo + vW-Vbo 



(22) 



From the above equation it is clear that for t so to be 
defined, 



M a + v (N - 3) b > 



(23) 



In other words, the continual collapse condition implies 
the positivity of the above term. Hence the time taken 
for other shells to reach the singularity can be given by 
the expansion, 



t s (r)=t S0 +rX{0) + O(r 2 ) 



(24) 



Also, from eauation (|17fl and l|19fl we get for small values 
of r, along constant v surfaces, 



v — v' = ^f(M + v( N -Vb ) (X(v) + rX 2 (v) + ■■■) 

(25) 

Now we can easily see that the value of ^(0) depends on 
and is completely characterized by the functions M(r) 
and b(r), which in turn specify fully the initial mass and 
energy distributions for the collapsing matter. Specifying 
these functions is equivalent to specifying the regular ini- 
tial data for collapse on the initial surface t = 0. In other 
words, a given set of density and energy distribution com- 
pletely determines the slope to the singularity curve at 
the origin, which is the central singularity. Also, it is 
evident that given any one of these two profiles we can 
always choose the other one in such a manner so that the 
quantity X(0) will be either positive or negative. 

In order to determine the visibility or otherwise of the 
central singularity, we now need to analyze the behaviour 
of non-spacelike curves in the vicinity of the singularity 
and the causal structure of the trapped surfaces. 

The boundary of the trapped surface region of the 
space-time is given by the apparent horizon within the 
collapsing cloud, which is given by the equation, 



F 



R N-3 



(26) 
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What we need to determine now is when there will be 
families of non-spacelike paths coming out of the singu- 
larity, reaching far away observers, and when there will 
be none. The visibility or other wise of the singularity 
is decided accordingly. Broadly, it can be stated that if 
the neighborhood of the center gets trapped earlier than 
the singularity, then it is covered, otherwise it is naked 
with families of non-spacelike future directed trajectories 
escaping away from it. By determining the nature of the 
singularity curve and its relation to the initial data, we 
are able to deduce whether the trapped surface formation 
in collapse takes place before or after the singularity. It is 
this causal structure that determines the possible emer- 
gence or otherwise of non-spacelike paths from the sin- 
gularity, and settles the final outcome in terms of either 
a BH or NS. 

To consider the possibility of existence of such families, 
and to examine the nature of the singularity occurring 
at R = 0, r = for the scenario under consideration, 
consider the outgoing radial null geodesies equation, 



dt 
dr 



(27) 



The singularity occurs at a point v(t s (r),r) = 0, which 
corresponds to R(t s (r),r) = 0. Therefore, if we have any 
future directed null geodesies terminating in the past at 
the singularity, we must have R — > as t — > t s . Now 
writing eauation (|27|l in terms of variables (it = r a ,R) 
where a > 1, we have, 



dR 
du 



-(«-!)# 



w 



(28) 



Choosing a = jf^j, and using eauation H13|) together with 
the collapse condition R < 0, we get, 




If null geodesies terminate at the singularity in the past 
with a definite tangent, then at the singularity we have 
4^ > 0, in the (u, R) plane which must have a finite 
value. 

In the case under consideration, all singularities for 

oo in the limit of approach 



r > are covered since ^ 



dR 



to the singularity in that case, and hence 
Therefore only the singularity at the central shell could 
be naked. 

In order to see possible emergence of null geodesies 
from the central singularity, we now need to analyze 
eauation (|29|l . Let us calculate the limits of the concerned 
functions in eauation H29fl at the central singularity. From 
equation® we get that in the limit of t — > t s ,r — > we 
have e~^R' — > 1. Also from eauation (|17(l and i|25|) . we 
have in this limit R — > 0. It then follows in general, from 
the Einstein equations discussed above, that the term 
F/R N ~ 3 goes to zero in this limit. 



We would like to find when there will be future di- 
rected null geodesies coming out from the central singu- 
larity with a well-defined and definite positive tangent 
in the (t, r) or (R, u) plane, thus making the singularity 
visible. Let us define the tangent to the null geodesic at 
the singularity as, 



,. ,. R dR 
.To = hm lim — = — — 

t-*t, r^a u du 



(30) 



*t,;r— >0 



Using equations l|29|l and l|25fl. along with the required 
limits as above, we get, 



(31) 



Let us now deduce the necessary and sufficient con- 
ditions for a naked singularity to exist, that is, for null 
geodesies with a well-defined tangent to come out from 
the central singularity. Suppose we have X(0) > 0, then 
we always have (from eauation (|31[l ). xo > and then in 
the (R, u) plane, the equation for the null geodesic that 
comes out from the singularity is given by 



R — xqu 



(32) 



In other words, equation is a solution of the null 
geodesic equation in the limit of the central singularity. 
Thus given X(0) > 0, we can always construct a solution 
of radially outgoing null geodesies emerging from the sin- 
gularity. This makes the central singularity visible. In 
the (t, r) plane, the null geodesies above near the singu- 
larity will be given as, 



t - t s (0) = x r~' 



(33) 



It follows that X(0) > implies xq > and we get 
radially outgoing null geodesies emerging from the sin- 
gularity, giving rise to the central naked singularity. 

On the other hand, if X(0) < 0, then we see that the 
singularity curve is a decreasing function of r. Hence the 
region around the center gets singular before the central 
shell, and after that it is no more in the spacetime. Now 
if there would have been any outgoing null geodesic from 
the central singularity, it must then go to a singular re- 
gion or outside the spacetime, which is impossible. Hence 
when ^(0) < 0, we always have a black hole solution. 

If X(0) = then we will have to take into account 
the next higher order non-zero term in the singularity 
curve equation, and do a similar analysis by choosing a 
different value of a in equation l|28|l . 

We have thus shown above that A'(O) > is the nec- 
essary and sufficient condition for null geodesies to come 
out from the central singularity with a definite positive 
tangent. It should be noted, however, that in general the 
dependence of R on r along the outgoing null geodesies 
from the singularity does not necessarily have to be of a 
power-law form. However, in order to satisfy the regu- 
larity and physical relevance, examining the trajectories 
which come out with a regular and well-defined tangent 
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is physically more appealing, which is the case we have 
examined here. 

We show below for completeness, however, that if null 
geodesies of any form come out at all, then those with 
definite tangent also must emerge from the central sin- 
gularity. Towards this, let us consider the equation of 
apparent horizon, which is from equations (|26|l and i|18|) 
given by, 



t ah (r) =t so + rX(0)+r 



..Mo) 



0{r 



iV-3 



(34) 



Since the apparent horizon is a well-behaved surface as 
one initiates close to r = 0, for a spherical dust collapse, 
hence we can say that the singularity curve for the col- 
lapse and the derivatives around the center are also well- 
defined, as the same coefficients are present in both (|34|) 
and (|19l) . Also, this shows that whenever ^(0) is neg- 
ative, the region around the center gets trapped before 
the central singularity, giving a sufficient condition for a 
black hole to develop. It follows that if null geodesies are 
coming out, then at least one of the coefficients X must 
be non- vanishing and positive. Then, as we have already 
shown, null geodesies with definite tangent will come out 
from the central singularity. 

From the above it follows that in the absence of a null 
geodesic with a definite tangent there cannot be any null 
geodesies coming out of the singularity. 

It is also clear now, from equation (20), that whether 
X(0) > or otherwise is fully determined by the regular 
initial data for collapse, in terms of the given initial den- 
sity and energy distribution for the collapsing shells. It 
thus follows that the initial data here completely deter- 
mines the final fate of collapse in terms of BH/NS end 
states. We shall discuss this further in the next section. 



IV. COSMIC CENSORSHIP AND GENERICITY 
ISSUES 

It is now possible to examine the question of valid- 
ity of cosmic censorship in a higher dimensional collapse 
scenario under consideration. Since we can explicitly find 
out as pointed out above, when families of non-spacelike 
geodesies can come out of the singularity or otherwise, 
we can address now this question below in some detail, 
and we also discuss some related issues. 

To focus the discussion, let us consider a model initial 
density profile as given by, 



14- \ I I 2^2 3 P3 

p{ti,r) = po + rpi + r — + r — 
and we write the function M(r) as, 

M{r) 



(35) 



n=0 



M n r r ' 



Mr 



2pn 



(N- 2)(N + n- l)n! 

(36) 

along with an energy profile as specified by, 



b(r) 



rbi+r — 



(37) 



A. Marginally bound collapse 

In the first place, let us consider the class of marginally 
bound collapse models for a more transparent under- 
standing of the problem. This is the case when the energy 
function b(r) above vanishes identically for the collapsing 
shells. In such a situation, the first non-vanishing coef- 
ficient X n (0), where n > 0, as described in equation 119|) 
are given by, 



X n (0) = - 



N -1 



Mn 



(38) 



The quantities M n are described in equation (|36|l . Now 
it is evident that whenever p± < 0, we will get a naked sin- 
gularity in all dimensions, whereas pi > always results 
in a black hole. The case p\ < corresponds to the phys- 
ical situation when the density decreases with increasing 
comoving radius r, as one would typically expect the den- 
sity to be highest at the center and then gradually de- 
crease as we move out in any realistic configuration such 
as a massive star. Further, we note that the above con- 
clusion is not dependent on the magnitude of p\ , but only 
on its sign, that is density must decrease away from the 
center with the density gradient being non-zero. Thus it 
becomes clear that it is the density inhomogeneity that 
delays the formation the trapped surfaces, thus causing 
a naked singularity. This is closely connected to the non- 
vanishing spacetime shear, and for a discussion of how 
shear distorts the geometry of trapped surfaces closer to 
the spacetime singularity, we refer to [l6| . 

Let us, however, assume that the initial density distri- 
bution has all odd terms in r vanishing, i.e. it admits 
no 'cusps' at the center and that it is either sufficiently 
differentiable, or is a smooth and analytic function of r. 
In that case, we must have p\ = 0. Then we get from 
equation ll'ol) that in the neighborhood of the singularity, 
the behavior of v is given by, 



lim lim v 



f N- 



(39) 



Also, in the same limit the function r n~3 has the form 



F 



r 2 M 



(40) 



Thus it is clear from equations l|39|) and (|40|l . that if 
N > 5 then for limt^t s , lim r ^o, ^ — ► o° an d thus the 
end state of collapse will always be a black hole (see also 
[10]). It thus follows that for a marginally bound dust 
collapse, with p\ — 0, i.e. when the initial density pro- 
file is sufficiently differentiable and smooth, the CCC is 
always respected in a higher dimensional spacetime with 
N = 6 or higher. 
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B. Critical value of p-2 in five-dimensional case 

We note that in the conclusion above, the spacetime 
dimension has to be either six, or higher. Let us consider 
the case when the spacetime dimension is five, but still 
with an analytic initial density profile. 

In case of a five-dimensional marginally bound collapse 
with pi = 0, we can write the tangent to the outgoing 
radial null geodesic at the singularity in the (R, u) plane 
as, 



1 - 



(41) 



The sufficient condition for the existence of an outgoing 
null geodesic from the singularity is we must have xq > 0, 
which in the above case amounts to, 



M 1< 



(42) 



But again, the outgoing null geodesic should be within 
the spacetime i.e the slope of the geodesic must be less 
than that of the singularity curve, 



lim lim ( — 

t->t„ r^o V dr 



null \ ' / sing 

From the above equation we get the condition, 

e- M *< 8 



(43) 



(44) 



Thus from equation (|42[1 and l|44|) we see that for an 
outgoing null geodesic from the singularity to exist we 
must have £ < £ c = —8, in which case we get a naked 
singularity, and otherwise a black hole results as collapse 
end state. 

We note that this situation has an interesting parallel 
to the four-dimensional collapse scenario, where we have 
a similar critical value existing, however, it is for the 
coefficient p3, when both pi and p2 are vanishing [l7| . 
Thus, with the increase of the spacetime dimension by 
one, the criticality separating the BH/NS phases shifts 
at the level of second density derivative from the earlier 
third density derivative. 

An interesting observation we could make here is that 
for £ < — 2 we have an increasing apparent horizon at 
the singularity. The apparent horizon is given by R = F 
so it initiates at the central singularity r = R = 0, and 
in the above case it is increasing in time (as opposed 
to the Oppenheimer-Snyder case of homogeneous dust 
collapse). Thus for the range — 2 > £ > — 8 no trapped 
surface is formed till the singularity epoch, but still we 
get a black hole as the collapse end state. This confirms 
that the absence of a trapped surface till the singularity is 
necessary, but not a sufficient condition for the formation 
of a naked singularity. 



C. Marginally bound collapse with initial 
homogeneous density 

It is useful to note here that the well-known 
Oppenheimer-Snyder class of collapse solutions is a spe- 
cial case of marginally bound dust collapse in four dimen- 
sions in which the initial density profile is homogeneous, 
that is, A4 n (n > 0) = for all n. The point is, if the 
initial density is homogeneous but if the collapse is not 
marginally bound, then the non-zero energy function / 
could inhomogenize the collapse at later epochs. In the 
present case, as / = 0, at all later epochs also the density 
remains a function of time only, that is, it is homogeneous 
at all later times as well. Thus we clearly see that the 
final outcome of this class of collapse is always a black 
hole. Furthermore from equations (|24|l and l|19|) we see 
that, 



t s {r) = t s 



,(t,r)=v(t) 



(45) 



As the scale function v is independent of r, all the shells 
collapse simultaneously to the singularity. The time 
taken to reach the singularity is given as, 



t. 



(N - 1)M§ 



(46) 



Thus, as we go to higher dimensions, for a given density, 
the time taken to reach the singularity will reduce. 

It is, however, interesting to note that even if we start 
with an initial homogeneous density profile but allow 
for non-zero initial radial and tangential pressures of the 
form, 



Pr(t,r) = 1 ;pe Q (r) = 1 + pe 2 r 2 +pe 3 r 3 + 
then as shown in [Tsj . we have 



X(0) 



(47) 



(48) 



Thus it is seen that a negative pe 3 coefficient does lead 
to a naked singularity. 

This is similar to the case of a collapse which is not 
marginally bound, where an initially homogeneous den- 
sity profile can turn inhomogeneous at later epochs due 
to the non-vanishing shell velocities. In the same way, 
in the case above also the non-vanishing pressures could 
inhomogenize the initially homogeneous density distribu- 
tion at later epochs to cause a naked singularity eventu- 
ally. It has to be noted, all the same, that the equation of 
state in the situations such as above could be considered 
to be some what peculiar (though matter is fully normal, 
satisfying the positivity of energy condition, and collapse 
conditions are fully regular). To state this differently, 
one can argue that the models where only purely tangen- 
tial pressures are taken to be non-vanishing may not be 
considered to be physically realistic, and if we choose the 
equation of state to be say p = hp, k > 0, or any homo- 
geneous equation of state, then when the initial density 
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profile is taken to be homogeneous, then so will be the 
initial pressures, and then the collapse will end up in a 
black hole only, and no naked singularity will arise. 



D. Genericity and stability of naked singularities 

In the above, we have discussed various special sub- 
cases of a higher dimensional collapse scenario which re- 
sult either in a black hole or a naked singularity end state, 
depending on values and behaviour of the parameters in- 
volved, and we also discussed the possible preservation 
of CCC under various set of assumptions. 

It is necessary however, to look at the situation in a 
collective manner if we are to gain any insight on the 
genericity and stability aspects connected to the naked 
singularities forming in gravitational collapse. It is well- 
known of course, that the genericity and stability are 
quite involved issues in the general theory of relativity, 
and that there does not exist any well-defined way to 
test the same in a unique manner for a given situation. 
Also, there may be different kinds of stabilities involved. 
For example, we can ask here, if the conclusions above 
will be stable to non-spherical perturbations, or when 
forms of matter more general than dust are considered 
and so on. Such issues are worth a separate and detailed 
investigation, and will be crucial towards the important 
problem of collapse end states. 

At a somewhat different, but still quite interesting 
level, we can inquire about the stability of BH/NS end 
states with respect to the perturbations in the initial data 
space which determines the final outcome of collapse as 
we have seen above. As pointed out here, this is a func- 
tion space consisting of all possible mass functions F and 
energy functions /. It is worth knowing how, for exam- 
ple, a naked singularity end state would be affected when 
one moves from a given density and energy profile (which 
gave rise to this state) to a nearby density or energy pro- 
file in this space of all initial data. 

The issue of how given density and energy distributions 
determine the final collapse end state has been discussed 
quite extensively in the usual four-dimensional dust col- 
lapse models, though using a somewhat different method- 
ology [l^. These results were completed to give a full 
and general treatment of four-dimensional case in [l9| . 
and the typical result is that given any density profile, 
one could choose the energy profile (and vice- versa), so 
that the collapse end state would be either a black hole, 
or a naked singularity depending on this choice. 

Our results here generalize this to the case of a higher 
dimensional collapse situation, and our method now al- 
lows us to make a more definite statement on genericity of 
naked singularity formation (in the sense stated above). 
As we see from equation (20), the quantity X(0) is fully 
determined from the initial data functions and their first 
derivatives. Once it is positive, the collapse ends in a 
naked singularity and a negative value gives black hole. 
It follows by continuity that given a density profile, if 



the energy profile chosen is such that the collapse ends 
in NS, i.e. X(0) > 0, then there is a whole family of near 
by velocities such that this will continue to be the case, 
and the NS will form an open subspace in the initial data 
space. Same of course holds for BH formations as well, 
and both these are neatly separated open regions in the 
initial data space. But if we take on physical grounds 
both p\ and b\ to be vanishing, then from equation (20) 
we have X(0) = 0, and CCC may be restored. 



V. CONCLUDING REMARKS 

In this final section, we give several concluding remarks 
and observations. 

1. It follows that if we can suitably motivate physically 
all the assumptions such as those discussed above, then 
it may be possible to restore CCC in a higher dimen- 
sional spacetime. Let us discuss such conditions in some 
detail. That the equation of state must be dust-like in 
the final phases of collapse is a strong assumption, but it 
is not a possibility that can be completely ruled out (see 
e.g. [13]). After all, we know very little on the equations 
of state, especially how it would be like in the advanced 
stages of collapse. Also, it is quite possible that in the 
very late stages of collapse the configuration is very much 
like a marginally bound one in the vicinity of the singu- 
larity. The introduction of pressures may or may not 
change such a scenario. 

All the same, in our view, the assumption p\ — is a 
tricky one, and this has been extensively discussed in the 
past within the context of collapse in four-dimensions. 
While it may be quite convenient to deal with smooth 
and analytic density profiles, especially when it comes to 
numerical models, it should not be forgotten that after 
all this is only an extra assumption, and that the basic 
equations of general relativity do not demand any such 
constraint. Neither it is clear astrophysically that the in- 
teriors of the stars must necessarily have analytic density 
or energy distribution. In certain equilibrium cases, the 
field equations imply that these have to be smooth, but 
this need not be true in general, and especially the dy- 
namically developing collapse situations could be quite 
different. 

We thus conclude that each of the above assumptions 
require further scrutiny and sufficient physical motivation 
so as to arrive at any definite conclusion on the status of 
CCC in a higher dimensional spacetime. However, this 
will be certainly worth the effort, given the importance 
of these issues in black hole physics. 

2. Let us consider the scenarios when some of the 
above assumptions break down. We have already noted 
earlier that when p\ is non-vanishing, then the collapse 
ends in a naked singularity in all dimensions, including 
N = 4. Again, our considerations above immediately 
imply that whenever spacetime is not marginally bound, 
the collapse always results into both the BH/NS phases 
as collapse end states, depending on the nature of ini- 



9 



tial data, irrespective of p\ being either zero or non-zero. 
That is, in a more generic non-marginally bound case, 
the condition p\ — does not save the CCC. However, if 
we take on physical grounds both b\ = pi = then from 
equation (20) X(0) — and CCC could be preserved. 

3. Consider the situation when we must believe some 
how that in the later stages of collapse the form of matter 
cannot be dust-like, and that non-dust forms of matter, 
and effects of pressures must be suitably taken into ac- 
count. In such ci CRSG, clS pointed out above, it is seen that 
even if one considered only homogeneous initial density 
profiles (however, with non-zero initial pressures), then 
also the pressure by it-self can cause sufficient distortions 
in the formation of the apparent horizon so as to cause 
a naked singularity as end state of collapse, rather than 
a black hole. However, if the equation of state is homo- 
geneous, together with initial data being homogeneous 
density profile, then no naked singularity will appear. 

4. We would like to suggest that there may be some 
hope, as outlined above, to recover CCC while we transit 
to a higher dimensional spacetime arena. This is subject 
to validity of several extra physical inputs as we described 
above. On the other hand, once we move to more gen- 
eral situations of either a non-marginally bound case, or 
with a more general form of matter, or without restric- 



tive extra-assumptions on the nature of the initial density 
profiles, then generically both the BH/NS phases could 
result as end states of collapse in a higher dimensional 
spacetime scenario. It would be fair to state that dynam- 
ical collapse in general relativity offers a rich spectrum 
of possibilities to investigate. 

5. Finally, we point out that the formalism given here 
brings out the role of initial data in causing BH/NS end 
states for the four-dimensional dust collapse also in a 
clear and transparent manner, completing earlier results 
in this direction. To be specific, it is seen, using equation 
(20), that given any initial density distribution for the 
cloud, one can always choose a suitable energy profile 
so that the evolution could end in either of a black hole 
or a naked singularity depending on the choice made. 
In other words, there is a non-zero measure of energy 
distributions which will take the given density profile to 
a black hole, and the same holds for a naked singularity to 
evolve from the same initial density. The converse is also 
true, namely, given any initial energy distribution, one 
can choose the density profiles which give rise to either 
of these end states. 
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